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Abstract 

The Lagrangian complex-space singularities of the steady Eulerian flow with 
stream function sin x\ cos X2 are studied by numerical and analytical meth- 
ods. The Lagrangian singular manifold is analytic. Its minimum distance 
from the real domain decreases logarithmically at short times and exponen- 
tially at large times. 
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1 Introduction 

Singularities of the solutions of hydrodynamical problems in the complex space 
domain are important, for example, because they allow us to give an objective 
definition of the otherwise fuzzy concept of "smallest scale present in a flow". 
Specifically, if a flow is analytic in the space variable x for all x, it can be con- 
tinued to complex z = x + iy locations where it will generally have singularities 
on some complex set S. The minimum distance 5 of £ to the real domain, called 
the width of the analyticity strip, defines the smallest scale. Indeed, the modu- 
lus of the Fourier transforms of the velocity v decreases roughly as e~~ 5k at high 
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wavenumbers k; thus the mesh needed in numerically simulating such a flow, e.g. 
by spectral methods, has to be significantly less than 5 (Sulem, Sulem and Frisch, 
1983; Brachet, Meiron, Orszag, Nickel, Morf and Frisch, 1983). Furthermore, 
should the flow ever develop a real singularity at a finite real time it must be 
preceded by complex singularities at a distance 5(t) which continuously vanishes 
att* (Bardos, Benachour and Zerner, 1976; Benachour, 1976a,b). For incompress- 
ible flow actual measurements of S(t), using high-resolution spectral simulations, 
indicate that 6(t) decreases in time, but in a way much tamer than suggested by 
the known rigorous estimates. In both two and three dimensions the behavior of 
8(t) at large times may well be exponential, but the best estimates are a decreasing 
double exponential in two dimensions and finite-time blow-up in three dimensions 
(see Frisch, Matsumoto and Bee (2003) for review). 

The standard explanation for this discrepancy is the phenomenon of nonlinear 
depletion, whereby the flow is found to organize itself into structures which have 
nearly vanishing nonlinearities (Frisch, 1995; Majda and Bertozzi, 2002). Note 
that one of the features characterizing depletion can be the degree of velocity- 
vorticity alignment or "Beltramization" of the flow. 

Up to now singularities in the complex domain were studied in Eulerian co- 
ordinates. In the present paper we analyse some aspects related to singular be- 
haviour of flows extended into the complex domain in Lagrangian coordinates. 
One practical motivation of studying Lagrangian singularities in the complex do- 
main concerns the structure of Eulerian singularities of a passive scalar. Indeed, 
consider a passive scalar 3 advected by an incompressible flow v, whose density 
satisfies the continuity equation 

d t 6(x,t) + V-(9v) = 0, (1) 

as long as we can neglect molecular diffusion. Obviously, this can be solved by 
using Lagrangian coordinates: 

0(x,t) = e o (a(x,t)), (2) 

where 9 is the initial density field and a(x,t) is the inverse of the Lagrangian 
map x (a, t), the latter satisfying the characteristic equation 

x = v(x, t), x(a, 0) = a. (3) 

If 6q is devoid of singularities (entire function), it is clear that the singularities 
of the Eulerian field, continued to complex locations, 9(z, t), will correspond to 

3 Similar remarks can be made about a passive magnetic field. 



2 



those (complex) fluid particles trajectories which have started at t = at (com- 
plex) infinity. Furthermore, it implies that the width of the analyticity strip 8(t) 
for 9(x, t) in Eulerian coordinates equals the width of the analyticity strip 5i,(t) 
of the Lagrangian map x(a, t) for the time reversed flow. In Section 4 we will use 
this correspondence to define the smallest scale of the passive scalar field 9(x, t) 
at time t. 

Note that the above argument can be applied to the 2-D vorticity. Assuming 
infiniteness of Eulerian vorticity at complex Eulerian singularities (Frisch, Mat- 
sumoto and Bee, 2003; Matsumoto, Bee and Frisch, 2003) it implies that the Eule- 
rian singularities of solutions of the 2-D Euler equation starting from entire initial 
data come from the complex infinity. Therefore, especially for the 2-D flows, 
some understanding of complex singularities in Lagrangian coordinates, even for 
a very simple flow, can shed light on the nature of Eulerian complex singularities. 

S. Orszag (2003) pointed out that the phenomenon of depletion is intrinsically 
Eulerian and has no Lagrangian counterpart 4 . Thus, it is expected that Lagrangian 
singularities are stronger - or at least closer to the real domain - than Eulerian 
ones. An extreme form of depletion 5 is to work with a Beltrami flow given by 
a steady solution of the incompressible Euler equation such as the ABC flow in 
three dimensions or the two-dimensional cellular flow with stream function 

if)(xi, X2) — sinxi COSX2 (4) 

and velocity 

V\ = — sin^i sin £2, ^2 = — cos X\ cos X2 . (5) 

The latter is much simpler, because the Lagrangian map can be obtained explicitly 
in terms of elliptic functions (Dombre, Frisch, Greene, Henon, Mehr and Soward, 
1986, Appendix A). Moreover, it represents a special case (A — 1, B — 1, and 
C — 0) of the ABC flow, obtained by applying a transformation of coordinates 
of the type used in Dombre, Frisch, Greene, Henon, Mehr and Soward (1986, 
Appendix A) and reducing the nontrivial dynamics to two dimensions. Our goal 
in this paper will be to obtain for this flow the complex singularities in Lagrangian 
coordinates. 

4 There are indications that depletion in the 3-D Lagrangian coordinates exists (Ohkitani, 2002), 
nevertheless it seems to be weaker than in the Eulerian coordinates. Furthermore, our definition of 
depletion (Frisch, 1995) is somewhat different from Ohkitani's definition of depletion, measuring 
the preference of the vorticity for being aligned with a certain eigendirection of the rate-of-strain 
tensor (the Eulerian case) and the Cauchy-Green tensor (the Lagrangian case). 

5 Observe that the nonlinear terms of the 2-D and 3-D Euler equations completely vanish for 
the flows (4) and the ABC flow, respectively. 
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The Lagrangian velocity 6 is defined just by the change of variables, that is, as 
v L (a, t) = v(x(a, t)) where x(a, t) is the solution of (3) with the velocity given 
by (5). Of course, such flow being steady, completely lacks Eulerian singularities 
(except at complex infinity). However, as we shall see, it develops real Lagrangian 
singularities in infinite time. 

The outline of the paper is as follows. In Section 2, we explicitly construct the 
singular set, the explicit construction of the Lagrangian map being relegated to 
Appendix A. In Section 3 the Lagrangian velocity, its Fourier transform and S(t) 
are determined numerically. Section 4 is devoted to the short-time and long-time 
asymptotics. 



2 Analytic structure of the singular manifold 

We consider the cellular flow on T 2 = [0, 2tt] x [0, 2tt] defined by (4) and (5). 
Because of the symmetries of this flow the full periodicity domain can be decom- 
posed into four cells, Fig. 1 (a). For convenience, in Fig. 1 (b), we represent the 
cell [0, 7r] x [— 7r/2, 7r/2]. The Lagrangian map x(a, t) is the solution of 

f Xi — — sin xi sin x 2 \ ,~ 
{ ■ (6) 

[ X2 = — COS X\ COS X2, 

with the initial conditions 

f xi(o, 0) = ai; 
\ x 2 (a, 0) = a 2 . 

When fluid particles start at complex locations they will stay in the complex do- 
main. We then denote the complexified Lagrangian position by c = a + ib and 
the complexified Eulerian position by z = x + iy. The complex Lagrangian 
stream function is £(c,t) = ip(z(c,t)) and the complex Lagrangian velocity is 
t> L (c,f) = v(z(c,t)). 

Since the Eulerian velocity v is an entire function of z , 7 the Lagrangian ve- 
locity can become singular only where the Lagrangian map becomes singular; 
furthermore it is easily shown that the only way z(c, t) can become singular, for 
a finite real t and a complex c, is by going to (complex) infinity. 

6 We cannot work with the Lagrangian vorticity since it remains constant along trajectories and 
thus constant in Lagrangian coordinates; the same applies to the stream function for any steady 
two-dimensional flow. 

7 An entire function is analytic for all z but can still be singular at infinity. 
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Figure 1 : Streamlines for the flow defined by (4) in the (a) flow box, (b) in the 
basic cell. 

The explicit solution to (6) is given in the Appendix in terms of elliptic func- 
tions. Its form implies the existence of singularities, associated to certain poles 
of suitable elliptic functions. This is however not the simplest way to actually 
construct the singular set a time-dependent complex manifold of complex 
dimension one, that is which can be parametrized in terms of one complex param- 
eter. For the construction it is simpler to observe (i) that the value of the complex 
stream function £ does not change along fluid particle trajectories and (ii) that the 
stream function, being the product of w\ = sin z\ and w 2 = cos z 2 , the only way z 
can run to infinity while conserving the stream function is to have either w\ — > oo 
and w 2 — > or conversely. 

In terms of the wi and w 2 variables, the equations for fluid particle trajectories 
can be rewritten as 




(8) 



with now complex initial conditions 



Wi(0) = sin c\ 



w 2 (0) = cosc 2 , 



(9) 



corresponding to the initial conditions (7). 



Obviously, the product £ = Wiw 2 is an integral of motion: this expresses just 
that for steady flow particles move along streamlines. One can take advantage of 
this integral of motion to rewrite (8) as a single equation in terms for example of 
the variable u = as 



u 



-Wl-u 2 ^l-£ 2 u 2 , u(0) = u = (1/Osinci. (10) 



This equation allows to express the time variable t in terms of the initial and 
current values of u by an elliptical integral: 

it = - \ -========. (id 



it = - 

J u 



y/(l-S*)(l-Ps*y 



We turn now to the determination of the singularities, that is those initial com- 
plex locations c* = (c*, c 2 ) which go to infinity at time t. Because of the symme- 
tries of the flow we can, without loss of generality, assume that when z — > oo at 
time t, we have w 2 — > oo and u>i — > and hence = 0. It follows then from 
(10) that the singular values of u are such that 

f° ds 
it = - . (12) 

J uo v/(l- s 2)(l-£V) 

The inverse of the elliptical integral is the elliptic sn function (Abramovitz and 
Stegun, 1965). Thus u can be expressed as an elliptic function. Since, by (10), 
wi(0) = £m , we obtain 

Wl (0) = £ sa(it, 0, w 2 (0) = ns(it, £), (13) 

where ns = 1/ sn. Alternatively, since the sn and ns functions of imaginary argu- 
ments can be expressed in terms of the sc function, defined as sc = sn/cn, and its 
reciprocal cs, we have 

w 1 (0)=itac(t,y/l-e), w 2 (0) = -ics(t, Vl-e 2 ). (14) 
By (9), we can express the singular Lagrangian locations at time t as 



c x = arcsm 



£sc(t, a/1 — £ 2 ) , c* 2 = arccos ics(t, a/1 — £ 2 ) 



(15) 



This constitutes the parametric representation of the singular manifold S(t), the 
parameter being the complex stream function £. Obviously, the functions appear- 
ing in this representation are multivalued. More precisely, the singular manifold 
has infinitely many sheets. 
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It is well known that elliptic functions are real for real values of the parameters. 
Therefore, S(t) will have a non-trivial part "above" the hyperbolic stagnation 
points, namely the four corners in Fig. 1(b). 8 By symmetry, it is enough to 
consider one of them. We take (ai, a 2 ) = (0, n/2) and easily obtain the following 
parametric representation of the singular locations having this stagnation point as 
real part: 



b\ = arcsinh £sc(i, a/1 — £ 2 ) 



arcsinh 



(16) 



Using the parametrization (16) in Fig. 2 (a) we have represented the curve 
obtained as intersection of the singular manifold with the set of points with real 
coordinates fixed at (x 1 ,x 2 ) = (0, tc/2) , that is the trace of the singular manifold 
on the pure imaginary plane lying above the corresponding stagnation point. Here 
"trace on" is used with its mathematical meaning of "intersection with". 
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Figure 2: (a) The cut of the singular manifold on the imaginary (y 1 , y 2 ) -plane 
is represented. The point closest to the origin lies on the y 2 -axis. (b) Here the 
value of the closest distance from singular manifold to the real domain, S(t), is 
represented for both positive and negative times. 

The width of the Lagrangian analyticity strip Si,(t) is obtained by finding the 
point(s) on the singular manifold closest to the real domain. It seems likely that 

8 By "above" we mean at complex locations whose real parts are the stagnation points. 
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these points will be located over the points where the flow exhibits a special struc- 
ture. Indeed, expanding (15) locally in its Taylor series confirms that the points 
closest to the real domain lie above the hyperbolic stagnation points. The upper 
left and lower right stagnation points have b\ = 0, while the upper right and lower 
left stagnation points have b\ = 0. The distance from the real domain is easily 
shown to be given by 



Fig. 2 (b) is a plot of S^t) showing both positive and negative times. 

As we have mentioned in the introduction, 5 L (t) of the time reversed flow 
gives the width of the analyticity strip in Eulerian coordinates for a passive scalar 
satisfying (1) with the same velocity field (5). Indeed, singularities for the passive 
scalar correspond to (complex) fluid particles which are mapped back to infinity 
in time t. Note that the symmetry of the cellular flow used here implies that S^t) 
is an even function of time so that S L (t) of the original flow coincides with the 
Eulerian 5(t) for the passive scalar field. 

3 Numerical integration in Lagrangian coordinates 

In this section we obtain the width 5 L (£) of the Lagrangian analyticity strip by 
numerically calculating the Fourier transform of the velocity in Lagrangian co- 
ordinates and then applying the method of tracing complex singularities (Sulem, 
Sulem and Frisch, 1983) to relate the Fourier transform to 5^(t). Our method 
is implemented for the steady cellular flow given by (4) but can in principle be 
applied to any steady solution of the Euler equation in both two and three dimen- 
sions. 

We need to calculate the Lagrangian Fourier coefficients, which are here ob- 
tained as follows. First, we calculate the Lagrangian map x(a, t) by solving the 
characteristic equations (6) with a fourth order Runge-Kutta method for N 2 ini- 
tial conditions a = ((2n/N)l, (27r/N)m) (I, m = 0, 1, . . . , N - 1), which form 
a regular square grid in the Lagrangian marker space. We obtain then the La- 
grangian velocity field by just changing variables: 



where v(x) is given by (5). The Lagrangian velocity Fourier coefficients are given 




(17) 



v L (a,t) = v(x(a, t), t), 



(18) 
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by 



Mk, t) = -r^^«L(a,t)e 



—ika 



(19) 



a 



For measuring the width 5^, it is convenient to define the shell-summed amplitude 



fc<|fc|<fc+i 

The Lagrangian width 5 L is then obtained by fitting the shell- summed amplitude 
to an exponential with an algebraic prefactor: 



We first discuss how to obtain the width 5 L (t) at short times when it is very 
large (since it is infinite at t — 0). Obtaining 5^(t) with a double-precision calcu- 
lation is not feasible because the Fourier amplitudes A^ L ^ fall off too quickly as a 
function of the wavenumber and thus get lost in the roundoff noise. So we need to 
solve the equation (6) with higher precision. 9 In practice we divide the time in- 
terval 1CT 12 < t < t into four parts: (i) 10~ 12 < t < 1CT 9 , (ii) 1CT 9 < t < 1CT 6 , 
(iii)lCT 6 < t < 10~ 3 , (iv) 10 -3 < t < 1. For each part, we use different time 
step At and precision: (i) At = 10~ 12 and 48-digit precision, (ii) At = 10~ 9 and 
36-digit precision, (iii) At = 1CT 6 and 24-digit precision, (iv) At = 10~ 3 and 15- 
digit (double) precision. These precisions, for each At, are the best we can handle 
with the fourth-order Runge-Kutta scheme. Note that in Frisch, Matsumoto and 
Bee (2003) a 90-digit spectral calculation was used for obtaining the short-time 
behavior of 5(t) in Eulerian coordinates for 2-D Euler flow with non-trivial Eule- 
rian dynamics; in the present case, the accuracy of the shell-summed amplitudes 
is constrained by the accuracy of the time integration scheme of (6). 10 Also 
we encounter the difficulty that a long-time integration causes serious accumula- 
tion of error in vorticity, which eventually breaks the conservation of vorticity (in 

9 We use here the package MPFUN90 (Bailey, 1995). 

10 We cannot take higher precision than the order of ( At) 4 , which is the error level of the fourth- 
order Runge-Kutta method employed here. If we take higher precision than that, we find that a 
bump appears in the shell-summed amplitude around the level of (At) 4 . This is perhaps due 
to the fact that the equation (6) is integrated in the physical space. In Frisch, Matsumoto and Bee 
(2003), the Euler equation was integrated with the same fourth order Runge-Kutta scheme but in 
the Fourier space, and this may considerably decrease the error (more precisely, the factor in front 
of (At) 4 is exponentially small for large k). 



A { ^> ofvL(Jfe) as 




(20) 



A^> oc AT a exp(-<S L £;). 



(21) 
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2-D vorticity should be conserved along each fluid particle trajectory). That is the 
reason why we split the time interval 1CT 12 < t < 1 into four parts and use differ- 
ent At for each part to avoid this accumulation. Figure 3 shows the shell-summed 
amplitude and the width of the Lagrangian analyticity strip 5i,(t) at short 
times 1CT 12 < t < 1. As seen in Fig. 3 (a), the number of points in to be 
fitted is small. We checked that, for several measured 5l(£)'s shown in Fig. 3 (b), 
a calculation with smaller At and higher precision for the same instance gives the 
same value of 5i(t). In this sense, we believe the measured is reliable in 
spite of the very limited data points. From Fig. 3 (b), a clean logarithmic decay of 
S L (t) is observed for more than 10 decades. 

Now we turn to the numerical determination of the long-time behavior of the 
width of the Lagrangian analyticity strip. The calculation is straightforward: the 
characteristic equation (6) is solved with standard double precision using a time 
step At = 10~ 3 . The number of grid points is 4096 2 . Figure 4 shows A^ L ^ and 
8 L (t) at long times, < t < 5. The temporal decrease of 5 L (t) is exponential in 
this regime. 

This indicates that exponentially small-scale structure is generated in the La- 
grangian coordinates by this steady flow. For seeing this, we first show in Fig. 5 
the deformation of the initially regular grid caused by the flow. As seen in Fig. 5, 
the Eulerian field is squeezed around the hyperbolic stagnation points along the 
contracting direction. We then plot contours of the modulus of the corresponding 
displacement \x(a, t) — a\ in the Lagrangian coordinates for various instances in 
Fig. 6. Here we can clearly see that the small-scale structures centered at hyper- 
bolic stagnation points are generated in the Lagrangian coordinates. The expo- 
nential decrease of S L (t) corresponds to the exponential decrease of the width of 
the structure in the Lagrangian coordinates centered at the hyperbolic stagnation 
points. 

4 Short-time and long-time behavior 

The parametric representation (15) for the singular manifold at arbitrary times is 
somewhat cumbersome. At short times we can use suitable expansions of the 
trigonometric and elliptic functions (Abramovitz and Stegun, 1965) and obtain to 
leading order 





(22) 
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Figure 3: (a) Shell- summed amplitudes of the Lagrangian velocity for four 
different short times calculated with the number of grid points 32 2 and 64 2 . (b) 
Log-linear plot of the width of the Lagrangian analyticity strip 5 L (t) measured as 
the logarithmic decrement of at short times 10 -12 < £ < 1. 
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Figure 4: (a) Shell- summed amplitudes of the Lagrangian velocity A ( ^ h> for 10 
different instances (t = 0.5, 1.0, . . . ,5.0) calculated with the number of grid 
points 4096 2 . (b) Linear- log plot of the width of the Lagrangian analyticity strip 
5 L (i) measured as the logarithmic decrement of at long times t < 5. 
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Figure 5: (a) Regular grids in the Eulerian coordinates (x 1 , x 2 ) at time t — 0. 
(b) Advected regular grids by the steady flow at time t = 1.6, which indicate the 
Lagrangian map x(a, t) in the Eulerian coordinates, (c) Advected regular grids at 
time t = 2.4. 
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Figure 6: (a) Contours of the modulus of the displacement \x(a, t) — a\ at time t = 
1.6 in the Lagrangian coordinates (a±, a 2 ). Contour levels are 7r/15, 27r/15, . . .. 
(b) Same as (a) but for t = 2.4. (c) Same as (a) but for t = 3.2. 
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from which it follows readily that the width of the analyticity strip is 

<5 L (0-lny, f->0. (23) 

Hence the Lagrangian 5 L (i) follows a logarithmic law at short time. 

We now observe that a short-time logarithmic law for the Eulerian S(t) has 
been obtained by Frisch, Matsumoto and Bee (2003) (see also Matsumoto, Bee 
and Frisch (2003)) for 2-D flow with a finite number of initial Fourier harmon- 
ics. In the present case the Eulerian dynamics are trivial since the flow is time- 
independent. We believe that a logarithmic law for the short-time behavior of the 
Lagrangian 5i,(t) is quite general for steady solutions of the Euler equation with a 
finite number of harmonics in any dimension d > 2. Roughly, the argument is as 
follows. Let n be the degree of the highest-order Fourier harmonic; the velocity 
at complex locations a large distance y from the real domain grows as e ny . Sin- 
gularities in Lagrangian coordinates at time t correspond to fluid particles which 
emanated initially from infinity or, equivalently, which are mapped to infinity in 
a time t under the reversed flow. The latter also grows as e ny at large y. A fluid 
particle located initially within a distance 5 L ^> 1 of the real domain and having 
an imaginary component of the velocity O (e ny ) will escape to infinity in a time 
approximately given by 



f 1 

t = I e~ ny dy = -e~ n5h . (24) 

It follows that 



5l 



n 



5 L (t)^-\n-. (25) 

n nt 

It is not clear whether this law for the Lagrangian 5 L will carry over to flow having 
non-trivial Eulerian dynamics. 

We now turn to the long-time behavior. In Eulerian coordinates for non-trivial 
dynamics the best proven lower bound for 5(t) is a decreasing double exponential 
(see Frisch, Matsumoto and Bee (2003) and references therein). But numerical 
simulations usually give a simple exponential decrease (Sulem, Sulem and Frisch, 
1983; Frisch, Matsumoto and Bee, 2003). The Lagrangian Si, for the steady cellu- 
lar flow considered here is given by (17), which has an obvious expansion at large 
t 

S L (t) - 2e~\ t -> oo, (26) 

which agrees with the numerical simulations of Section 3. Note that the nearest 
singularities are located "above" the hyperbolic stagnation points of the flow (5). 
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As we have remarked in Section 1 S L (t) allows us to give an objective defini- 
tion of the smallest scale for the passive scalar advected by the flow (5). According 
to formula (26) this scale decreases exponentially in time which, as we have noted 
in Section 3, reflects the squeezing of the flow around the hyperbolic stagnation 
points. The exponential temporal rate of deformation of passive scalar structures 
transversal to the streamlines can be explained by dimensional reasoning. 

We note that the period of motion of a particle along a trajectory with the 
stream function ^(a) is given by 4K(^1 — \l/ 2 (a)) (see Appendix A, formula 
(A.4) and (Lawden, 1989)) changing from 2n in the centre of the cell to infinity 
at the boundaries. Obviously, for a passive scalar field which is not constant 11 
along the streamlines, such as for example 9(x 1 ,x 2 ) = cosxi sinx 2 , the most 
intensive mixing will be observed near the hyperbolic stagnation points. Let us 
consider the flow near the stagnation point (0, 7r/2) atthepointa = (ai, ir/2— a 2 ). 
The displacement of a fluid particle starting at a becomes significant at the time 
t ~ K(y/l - ^ 2 (a)). Since *(o) is small, we use the approximative expression 
(see (Whittaker and Watson, 1927)) 

#(^1 -*>))« In (27) 

and obtain 

a x a 2 ~ e~*. (28) 

Hence the smallest spatial scale for the passive scalar will decrease exponentially 
in time. 

In this paper, it is found that hyperbolic stagnation points can be a key in the 
analysis the width of the Lagrangian analyticity strip of the steady solution (4) to 
the 2-D Euler equation. Currently we are investigating the Lagrangian 5^(t) of 
the 3-D steady solutions (the ABC flows and more general Beltrami flows (Ma- 
jda and Bertozzi, 2002)) of the 3-D Euler equations. Since such flows can be 
chaotic, an interesting question is which stretching — caused by the chaos or the 
hyperbolic stagnation points — dominates. A preliminary result indicates that 
the nearest singularities are again above the hyperbolic stagnation points at large 
times. The details of the study of the 3-D steady flows will be reported elsewhere. 
It is also of interest to point out that for kinematic dynamos in three-dimensional 12 
steady flows, the fastest growth of the magnetic field is frequently observed near 

n If a passive scalar field is constant along the streamlines there will be no mixing. 
12 Because of the 2-D antidynamo theorem three dimensions are required. 
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hyperbolic stagnation points (Soward, 1994) despite the presence of the chaotic 
stretching. 

A question for further study remains whether the behaviour of the Lagrangian 
8h(t) does exhibit some structural stability with respect to hyperbolic stagnation 
points. 
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Appendix 

A Explicit construction of the Lagrangian map 

The system of ordinary differential equations (6) can be solved by a simple adap- 
tation of what is done for an integrable case of the ABC flow in (Dombre, Frisch, 
Greene, Henon, Mehr and Soward, 1986, Appendix A). For completeness we give 
the full derivation, which is quite elementary. 

Taking the derivative of (6) with respect to time we obtain a system of two 
decoupled differential equations 

2xi = sin 2xi , . 

2x2 = — sin 2x2, 
satisfying the initial conditions (7) and 

ii(a, 0) = — sinctisina2 
X2(a, 0) = — cosaicosa2- 

Eq. (A.l) is obviously the same as a set of pendulum equations describing non- 
linear oscillations around the origin for the rr 2 -variable and around n/2 for the 
rri-variable. Let us just consider the rr 2 -variable. The equation has the first inte- 
gral 

-Xo — - cos 2x 2 = - — - sin 2 a\ cos 2 a 2 = - ( 1 — 2^ 2 (a)) , (A. 3) 

2 2 4 42 4 v ' 
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which expresses the conservation of energy. From this equation we obtain by 
standard integration (Lawden, 1989) 



*+/ Tfz T2 , „ — Vl-^ 2 (a) 
Jo a/(1 - ^ 2 (a)) - sin 2 a; / 

(A.4) 

Obviously, since the poles of the sn function do not lie on the real axis (Lawden, 
1989; Abramovitz and Stegun, 1965), solutions of (6) are periodic and nonsingular 
for real initial conditions. The only possibility for an orbit to come across a pole 
and thereby to go to infinity is to start at a suitable complex location. 
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